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The effect of a periodic pinning array on the vortex state in a 2D superconductor at low tempera- 
tures is studied within the framework of the Ginzburg-Landau approach. It is shown that attractive 
interaction of vortex cores to a commensurate pin lattice stabilizes vortex solid phases with long 
range positional order against violent shear fluctuations. Exploiting a simple analytical method, 
based on the Landau orbitals description, we derive a rather detailed picture of the low tempera- 
tures vortex state phase diagram. It is predicted that for sufficiently clean samples application of 
an artificial periodic pinning array would enable one to directly detect the intrinsic shear stiffness 
anisotropy characterizing the ideal vortex lattice. 
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I. INTRODUCTION 

The nature of the vortex lattice melting transition in 
2D superconductors has been debated in the literature 
for many yearsi Early proposals^, based on the simi- 
larity to the Kosterlitz-Thoulcss-Halperin-Nelson- Young 
theory of melting in 2D solids^, have led to the conclu- 
sion that the melting transition is continuous. A weak 
first order melting transition was predicted more recently, 
however, by several Monte Carlo simulations^ using 
the Ginzburg-Landau (GL) theory. It has been shown 
recentlj*ii£i£ that shear motions of Bragg chains along 
the principal crystallographic axis of the vortex lattice 
cost a very small fraction of the SC condensation energy 
and are responsible for the low temperature vortex lattice 
melting. 

This intrinsic anisotropy of the vortex lattice with re- 
spect to shear stress can not be easily detected experi- 
mentally since the orientation of the principal axis with 
respect to the laboratory frame depends on the local pin- 
ning potential, which in real superconductors is usually 
produced by random distribution of pinning centers. In- 
direct experimental detection of this hidden anisotropy 
may be achieved by means of the small angle neutron 
scattering (SANS) technique, due to the ID nature of the 
effective thermal fluctuations in the vortex liquid state 
just above the melting point (see Ref*-). A direct detec- 
tion of this anisotropy (e.g. by means of SANS) could 
be possible if vortex solid phases with long range posi- 
tional order were stabilized against the random influence 
of pinning impurities. This can be achieved by expos- 
ing the SC sample to an artificial periodic pinning array 
and tuning the magnetic flux density to an integer mul- 
tiple of the pinning centers density. As will be shown 
in this paper, under certain conditions the artificial peri- 
odic pinning potential can stabilize weakly pinned vortex 
solid phases with long range positional order, which may 
exhibit the shear stiffness anisotropy characterizing the 



ideal vortex lattice. 

Vortex matter interacting with periodic pinning arrays 
is currently a subject of intense experimental - — and 
theoreticalii - — investigations. Developments of nano- 
cnginccring techniques, such as e-beam lithography, make 
it possible to fabricate well defined periodic arrays of sub- 
micron antidotes, or magnetic dots, in SC films with low 
intrinsic pinning, enabling to study the effect of well con- 
trolled artificial pinning centers. These experiments have 
shown that under certain conditions the underlying arti- 
ficial pinning centers can attract vortices very strongly, 
thus stabilizing vortex patterns with global translational 
symmetry against the randomize influence of the natural 
pinning centers. 

From theoretical point of view the utilization of an 
external periodic pinning potential provides a convenient 
tool for testing different models of the vortex state by 
simplifying considerably the model calculations. At the 
same time, however, the interplay between the vortex- 
vortex interactions, which favor hexagonal vortex lattice 
symmetry, and the underlying periodic potential can lead 
to a variety of vortex configurations, depending on the 
pinning strength, in which vortices detach from pinning 
centers to form more closely packed vortex patterns. 

As the interaction with a periodic substrate stabilizes 
the vortex system versus thermal fluctuations, it gener- 
ally increases the melting temperature. However, as we 
shall see in this paper, deviation from the ideal hexagonal 
symmetry due to pinning reduces the phase dependent 
interaction between vortex chains^, making them less en- 
during under thermal fluctuations. In the weak pinning 
limit, where depinned floating state can occur, the cor- 
responding phase diagram becomes rather complicated, 
due to the possibility of transitions between floating solid 
and pin solid phases^. 

In the present paper we study the influence of a peri- 
odic pinning substrate on the vortex state in 2D , extreme 
type II superconductors, at perpendicular high magnetic 
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fields. Our approach is based on the previously devel- 
oped theory of vortex lattice melting in pure supercon- 
ductors carried out within the framework of the 
GL theory in the lowest Landau level (LLL) approxima- 
tion. Specializing the calculation for a vortex system 
interacting with a square pinning array under the first 
matching magnetic field, we study in detail some key lim- 
iting regions of the vortex phase diagram, which enables 
us to determine its main qualitative features. 




II. THE MODEL 

We consider a 2D superconductor at high perpendicu- 
lar magnetic field, interacting with a periodic substrate 
of pinning centers, located at (xi,yj). A phenomeno- 
logical Ginzburg-Landau functional, with an order pa- 
rameter tjj (x, y), is used to describe the superconducting 
(SC) part of the free energy, and a local periodic pinning 
potential^: 

Vpin = v ^2 HixiiVj) I 2 (!) 

determines the interaction of the vortex state with pin- 
ning centers. We assume vq > 0, so that the pinning en- 
ergy is minimal if the vortex core positions, determined 
by ip (xi, yj) — 0, coincide with pinning centers. 

Our main interest here is in the influence of the pinning 
potential on the vortex lattice melting process, so that 
the pinning energy V P i„ is restricted to the range of the 
vortex lattice melting energy, which is much smaller than 
the SC condensation energy. Since the latter is of the 
same magnitude as the cyclotron energy, it is justified 
to restrict the analysis to the LLL of the corresponding 
SC order parameter, which can be therefore written as a 
linear combination of ground Landau orbitals: 

V> (x, y) = ^2 c n<t>q n (x, y) (2) 

n 

c n = \c n \e^; 4> q (x,y) = e 2i « x -^+tf 

where q n = qn, q = n/a x , and the amplitudes c n in 
the mean field approximation are related to the (spatial) 
mean square SC order parameter, Aq , through: \c n \ 2 = 

Cq = \J 2|^Aq. In our notations all space variables are 
measured in units of magnetic length. 

In this model, due to the Gaussian attenuation along 
the y-axis over a characteristic distance of the order of 
the magnetic length, the vortex cores (located at the ze- 
ros of ip (x, y) ) form a network of linear chains along 
the x-axis, each of which is determined mainly by a su- 
perposition of two neighboring Landau orbitals^. The 
parameter a x is therefore equal to the inter-vortex dis- 
tance within a chain, while n/a x is the inter-chain spac- 
ing in the y-direction (see Fig. It should be noted 
that deviations of Landau orbital (LO) amplitudes from 



FIG. 1: Schematic arrangement of a vortex lattice relative to 
the pin lattice. 



their mean field value, Co, resulting in strong local dis- 
tortions of the superfluid density and large increase of 
the corresponding free energy density, are neglected in 
comparison with variations of the phase variables, ip n , 
provided the orbital direction is selected to be along the 
principal crystallographic axisi^. 

We select the pinning centers to form a rectangular 
lattice 

(xi, yj) = (IJ + x ,l y j + y Q ) (3) 

where i,j = 0, ±1, The parameters Xq and yo deter- 
mine the relative position of the pin and vortex lattices. 
The nature of the vortex state in the presence of the 
pinning potential depends crucially on the ratio of the 
number of vortices, N = s/N x y/N, to the number of 
pinning centers, N p = N PtX x N py . Since the density of 
vortices depends on the external magnetic field strength 
H , one can tune this ratio by varying H. Of special in- 
terest are the matching fields H = H u , v = 1, 2, when 
the ratio n§ = N/N p = v is an integer. 

In matching fields one may distinguish between two 
different situations, when vortices are bound or unbound 
to pinning centers. If the pin lattice and the vortex lattice 
unit cells are commensurate along both x and y directions 
, i.e. l x = c x a x , l y = c v n/a x , with c x and c y being 
integers, the pinning energy is equal to zero, since all 
the vortices coincide with pinning centers. In all other 
cases of matching fields , c^Cj, = n x integer , non of 
the numbers c x > 1 and c y > 1 can be integer , and the 
lattice constants are incommensurate in both directions. 
It will be shown below that such a vortex configuration is 
in a floating state with respect to the pin lattice, similar 
to vortex states in mismatching magnetic fields. 

Using the LO representation, Eq.JSJ), of the SC order 
parameter in Eq. Q for the pinning energy, one may take 
advantage of the localized nature of the LOs and expand 

_ 2 

V P i n in the small parameter A = e 9 , which reflects 
the small overlap integral between adjacent orbitals con- 
tributing to the local superfluid density at the pinning 
centers. Retaining only dominant terms in A Eq. Q is 
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reduced to the form: 




. in 



vortex is, up to small terms of the order ~ e 27r , given 
by: 



V m 



N 



+ 2e~ q /2 u k+1/2 <S> kA 



V~J= ^ [<Il - »2 COS - tpk-l)] 



u fc+m/2 



iV, 



^ e -2( W + 9 (fc+m/2)) J 



— ^ cos (9?fe +m - tfk + 2mqxi) 



N„ x 



(4) 



where Vb = VoA^Np^Np^y. It should be stressed that this 
approximation is valid only for LOs along the principal 
axes since the minimal distance between them, q = tt /a x , 
is sufficiently large to ensure small and rapidly decreas- 
ing value of the overlap integrals between more distant 
orbitals. 

If c x is not an integer, namely the rectangular pin lat- 
tice and vortex lattice are incommensurate in x-direction, 
then Eq.0} shows that &k,m = 0. In this case the pinning 
energy does not depend on the phases (i.e. the relative 
horizontal positions ) of the Landau orbitals. 

Expressing the functions Uk and Ufc+1/2 with the help 
of Poisson summation formula as 



Mfe+m/2 — 



— £< 

jv p,y a 



-2(l v j+y +q(k+.5)) 2 




we note that when the lattices are incommensurate also 
along the y axis ( i.e. when both c x and c y are not 
integer) the oscillating terms in u k are averaged to zero 
after summation over k. Thus, the pinning energy for 
incommensurate lattices is a constant 



V, 



pin 



(0) 



which does not depend on the mutual orientation of the 
vortex and the pin lattices. Note that the system size 
in y direction is L y = q\/~N = lyN PtV , a relation con- 
necting y/N to N P y. Obtained result is valid only for 
large system, N — ► 00, where the boundary effect can 
be neglected. 

For the sake of simplicity, we will consider in what 
follows a square pin lattice with 710 = 1. In the commen- 
surate situation the pinning energy is minimal ( i.e. equal 
to zero) when all vortices coincide with pinning centers. 
Deviations of vortices from this configuration in the form 
of shear distortions along the principal crystallographic 
axes are of special interest due to the relatively low SC 
energy involved. For the principal axis parallel to a side 



of the square unit cell, c x 



1 and q = tt, and 



where v = V /N , ax = 1 - 2e _7r / 2 ~ .584 and a 2 = 
2e~ 7T / 2 (l + 2e~ 7r / 2 ) ~ .589 . Note that in the above 
expressions we set xq — yo — q/2 so that the minimal 
pinning energy is obtained for = 0. Note also that for 
the undistorted square lattice in which tp^ = tpk—i, the 
expression in the first line of Eq. Q is not strictly zero 
since a\ ^ 02. The error, which is of order higher than 
the second in e~" , can be neglected in the approxima- 
tion leading to Eq. (JJJ). The numbers 01, 0.2 can be thus 
considered equal within this approximation, allowing us 
to introduce a single coefficient n x = a\ ~ ai ~ .59. The 
expression in the second line of Eq. Q yields the correct 
(i.e. zero) value for the undistorted lattice. It is writ- 
ten in terms of the variables, £fc = ipk — ifk-i, describing 
the lateral positions of the vortex chains, which are gen- 
erated mainly by interference between two neighboring 
LOs. This is consistent with the well known definition 
u x = dtp/dy of vortex displacement along the x axis in 
the long wavelength limi^. 

To evaluate the excess pinning energy associated with 
shear distortion along the diagonal of the square unit 
cell the pin lattice may be conveniently described by 
two interpenetrating simple square sub-lattices with c x = 
l,c y = 2 and q 2 = 2tt ( see Fig|2). The corresponding 
interchain pinning energy for each of the sublattices can 
again be obtained from Eq. l[7|l. with k x i ~ .84 and a 
phase shift of tt/2 , which arises due to different shape of 
the unit cell. 

The SC part of the free energy functional for the com- 
mensurate lattices described above (a x — y/n) is given 
by the following (^-dependent) expression^: 



H sc 



= -kr 



(8) 



where hn is the SC condensation energy (per unit flux) 

4A 2 

of the square vortex lattice, and T n = hn is the 

shear distortion energy parameter, expressed through 
the dimensionless interchain coupling constant, X sq = 



Here h L 



£oj^- , where /3a — 1.159 and 



so, according to Eq. the pinning energy per single 



CXp^ — 71 

(3 sq ~ 1.18 , are the values of the Abrikosov structure 
parameter for regular hexagonal and square lattices re- 
spectively, and £0 is the SC condensation energy of the 
former. 

For the specific choice = jk, where 7 is a constant, 
the Bragg family of vortex chains along the principal 
axis, denoted x, is characterized by a lateral displacement, 
Cfe+i ~ £fc = 7? between neighboring chains. Evidently, 
the SC energy, H sc , for the undistorted square vortex 
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FIG. 2: Primitive and none-primitive unit cell representations 
(solid and dotted lines respectively) of the square pin lattice 
used for description of shear distortion along the principal 
axes of the vortex lattice. 
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lattice £ k = ( 7 = ) ( see Eq. © ) , is equal to Nh n ■ 
However, the minimum of the SC energy with respect to 
the collective tilt angle parameter 7 is reached for a tri- 
angular vortex lattice, determined by = nk (7 = n), 
whose unit cell is an isosceles triangle with a base ( along 
x-axis) and a height equal to y/ir (see Fig. |2J. The cor- 
responding SC energy is equal to H sc /N = —hn — 2Tb . 
This value is lower than the SC energy of the square vor- 
tex lattice, and only slightly higher ( i.e. by ~ .45%) than 
the SC energy of the equilateral triangular (Abrikosov) 
lattice, H/\/N = — e . 



III. VORTEX STATES FOR THE LOWEST 
MATCHING FIELD 

A. Commensurate and incommensurate Ground 
states 

The competition between the pinning energy, Eq. J7J), 
which favors vortices approaching the pinning points on 
a square lattice, and the SC energy, Eq. ©, prefer- 
ring triangular lattice configuration, leads to 'frustrated' 
vortex structures, which depend on the relative pinning 
strength. 

At zero temperature they can be obtained by mini- 
mizing the total energy, consisting of the SC and pinning 
parts. Since in the LO representation each orbital is v^V- 
fold degenerate, the effective Hamiltonian is written as: 



/□ = 



H sc H~ ^pm 

TV 



Ta 



+ 4p[l- cos(&)] - [1 - cos(a+i - &)] } (9) 

where the parameter p = K x v n /4T n determines the 
strength of the pinning potential relative to the inter vor- 



o- 



O x C x 



-0- 



(b) 



FIG. 3: (a) The vortex lattice state with the lowest energy, 
which is commensurate with a square pinning lattice in the 
limit of zero pinning strength, (b) An alternative vortex lat- 
tice state, which may be favorable under square framework 
boundary conditions (see text). 



tex chain coupling. Under the constraints imposed by 
the requirement of commensurability between the vortex 
configuration and the pin lattice, the vortex chains are 
restricted to move laterally along the common x-axis of 
the underlying lattices (see Fig. The corresponding 
displacements, , may be separated into two groups, 
corresponding to even and odd vortex chains, as follows: 



61 for k = 21 
for k = 2l-l 



(10) 



so that 



fn = -^E{-2/i n /T n +4p[2-cos(^)-cos(0)] 

- [2 - cos(0; - Ci) - cos(0,_i - CO] } (11) 

The calculation may be greatly simplified if we assume 
that the stationary point values within each group are all 
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equal, that is: 9i = 9 C and Q — £ c . This restriction may 
be justified in the weak pinning regime < p < 1 , where 
the dominant SC energy part favors periodic triangular 
vortex structures, as shown in Fig. 

Substituting these values to Eq. Hll|) and minimizing 
the resulting functional one finds that the total energy 
has a minimum: 

mm (/□) = --!= [ 2 W - Pf + M (12) 

at 9 C = — Cc = $o where 

cos(9 )= P . (13) 

Thus, at zero pinning strength, p = 0, the ground state 
energy per unit flux, /□ = — 2Tb — hn, corresponds to a 
triangular vortex lattice configuration, 9 C = — Cc = tt/2 
, whereas in the opposite extreme, when p — 1 , the 
ground state energy, /□ = hn, corresponds to a square 
vortex lattice, 9 C = — £ c = , which coincides with the 
underlying pin lattice. It should be stressed, however, 
that due to the constraints imposed by the requirement 
of commensurability with the pin square lattice, the tri- 
angular vortex structure obtained in the zero pinning 
strength limit, is not the equilateral (Abrikosov) lattice 
(see FigOHi). This discontinuity indicates that the tran- 
sition to the depinned (floating) vortex lattice should be 
of the first order (see below). 

An illustration of the weak pinning ground state con- 
figuration is shown in Fig. [3^. It is seen that odd and 
even vortex chains are shifted in opposite directions sym- 
metrically with respect to the underlying substrate. The 
relative positions of the two lattices are determined by 
the strength of the pinning potential, Vo- In the zero 
pinning limit, Vq — ► 0, the vortices in odd (even) chains 
approach lattice points which are shifted laterally by a 
quarter of a lattice constant, -A x ( l x = ), in the 
positive (negative) sense with respect to the square pin 
lattice, forming isosceles triangular lattice. Note that 
the asymmetric configuration, shown in Fig|3j3, in which 
half of the vortex chains remain pinned to the underlying 
substrate, has energy — 2Tb (1 — 2p) — ha, which is only 
slightly (i.e. by a small, second order correction in p ) 
higher than the energy given by Eq. (|12|1 . Such an asym- 
metric configuration may become energetically favorable 
(see Refi2£iSi) due to, e.g. boundary conditions which 
are incompatible with the even-odd chain symmetry de- 
scribed above. 

At sufficiently weak pinning, when the pinning energy 
becomes comparable to the difference between the SC 
energies of the commensurate isosceles triangular vortex 
lattice with a 2 , = it , and the incommensurate equilateral 
triangular lattice with a 2 , = v3tt/2 (Abrikosov lattice), 
the latter is preferable. To show this note that the energy, 
— ft, a j of the equilateral triangular vortex lattice in the 
presence of incommensurate pin lattice is influenced only 
by the average pinning potential v, so that: 

4 

- h-A = so +v= -s + —pT a (14) 

fox 



Comparing this value with that obtained in Eq. (|12f> for 
the commensurate, isosceles triangular lattice, —2Tb (1 — 
p) 2 — /in, we find that for p < p c — .25 the floating equi- 
lateral triangular lattice is the lowest energy state. This 
critical point can be thus identified as a transition point 
from pinned (commensurate) solid to a floating (incom- 
mensurate) solid state. 

A second critical point exists in the strong pinning 
regime, i.e. at p = 1, as indicated by Eq. l|13|) . which 
has no real solution at any p > 1 . At this critical point 
the vortex lattice coincides with the square pin lattice ( 
i.e. 9q = in Eq. l|13fl ) and the pinning energy reaches 
its absolute minimum value (i.e. zero ). Since any fur- 
ther increase of the pinning strength above the critical 
value, p = 1, can not be compensated by the SC energy 
terms in Ea. (|llfl . the vortex configuration remains fixed 
at the square lattice structure for any p > 1. Thus, with 
increasing values of the parameter p , the ground state 
vortex configuration changes continuously from a trian- 
gular lattice at p ~ p c , into a square lattice at p = 1 , 
which does not changes with further increase of the pin- 
ning strength. This continuous transformation from a 
triangular lattice to a square lattice can be classified as 
a second order phase transition at p = 1 . 



B. Commensurate equilibrium states at finite 
temperature 

In the ideal vortex state at finite temperature thermal 
fluctuations associated with the low-lying shear excita- 
tions along the principal crystallographic axis destroy the 
long range phase coherence of the vortex state and lead 
to melting of the ideal vortex lattice at a temperature, 
T m , well below the mean field T c . This feature indicates 
an intrinsic anisotropy of the ideal vortex crystal^: The 
characteristic excitation energy for sliding vortex chains 
along the principal axis ( denoted by x ) parallel to a side 
of the unit cell is two orders of magnitude smaller than 
the SC condensation energy, and one order of magnitude 
smaller for fluctuations along the diagonal axis ( denoted 
by x'). For all other crystallographic orientations the 
shear energy is of the order of the SC condensation en- 
ergy. 

The nucleation of a SC crystallite can be established 
in such an ideal model by selecting boundary conditions 
which fix the position of a single vortex chain with re- 
spect to the laboratory frame. As shear fluctuations of 
parallel vortex chains diverge with the distance from the 
fixed chaini, a SC domain is restricted to nucleate only 
around a pinned chain, its transverse size shrinking to 
that of a single magnetic length as the temperature rises 
toward T m . For the sake of simplicity, we avoid here 
the complication associated with the discontinuous na- 
ture of the vortex lattice melting process, which involves 
two principal families of easily sliding Bragg chainsi, and 
restrict the analysis to a single family of vortex chains, 
i.e. that with the lowest crossover temperature T crr ^. A 
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FIG. 4: Pair correlation function between nearest Landau 
orbitals: (1) at strong pinning, as a function of the dimen- 
sionless inverse temperature r = ApTu/T, and (2) in the tri- 
angular Abrikosov vortex lattice with a single pinned chain , 
n = 0, as a function of f = 



meaningful definition of T cm (a x ) may invoke the phase 
correlation function, 



Cn',?i — ( 



(15) 



between Landau orbitals, n' and n , located near the 
fixed chain n = . Thus, melting of the entire vortex 
lattice occurs essentially when phase correlation between 
the nearest neighboring chains ( i.e. n = 1, n' = 2 ) 
closest to the fixed chain is significantly suppressed (e.g. 
by a factor of 1/2). In the p — ► limit we use the 
expression derived in Ref.(£) to find: 



C n >- 



(r) 



h(r) 
I (r)J V W) 



'i/2(rj 



Ta/T . Here the characteristic temperature, 
£o with A = cxp(— \^3tt/2), corresponds to 



4\ z 



1+4A 



with r 

T A ~ 

interaction between the principal LOs in the equilateral 
triangular lattice state. Note that the crossover between 
the vortex solid state at zero temperature, where f — ► oo 
, and C n ' = 2,n=i (t^ — >• 1 , and the vortex liquid state at 
high temperature, where r — > , and C„'=2,n=i (j) ~ * 
, occurs at about r ~ 1.5 , so that T cm (a x ) ~ .677a 
(see Fig. 0J. This crossover temperature is close to , 
though somewhat lower than the melting temperature, 
T m ~ 1.2T A ~ 2.8Tb, predicted in Ref.(i). 

The presence of the periodic pinning potential stabi- 
lizes the vortex lattice against the violent phase fluctua- 
tions discussed above. This effect is nicely demonstrated 
by the phase correlation function C n > in ( Ea. (|T3|l ), which 
controls the mean superfluid density (see Eq.(J5J) near 
the melting point. Assuming strong pinning, p> 1, and 
neglecting the small GL inter-vortex-chain coupling, the 
correlation function can be determined from the exprcs- 



Uk Jo d^ k e z( -fr,'-'p n ) e -'iprcos(i k ) 



cos(5 fc ) 



(16) 



Using the identity <p n = ^22=n Cfc) where the value of 
no can be found from boundary conditions which influ- 
ence only the global phase of the SC order parameter, we 
find that 



|An| 



/i(4pr) V 

exp(-|An|/8pr) for Apr > 1, (17) 



This result contrasts with the cor- 
relation function obtained in the pure stated, which has 
the asymptotic form 



C„/,„ cx exp ( -7jHAn| 2 ) for r > 1 



where n = ^- + =j- — ^ , with the n = chain being 
fixed. As discussed above (see also Refs.(— *i^)), fixing 
chain positions through boundary conditions is physically 
equivalent to introducing pinning potential into the GL 
functional, which is a crucial step for stabilizing the vor- 
tex lattice. The global stability of the vortex lattice in the 
presence of the periodic pinning potential is reflected in 
Ea. (|17|) . as compared with Eq.JTSJ, by the translational 
invariance of the former correlation function, as well as 
by its relatively weak (simple exponential) decay. 

To determine the crossover temperature from the 
square pin solid (SPS) to the vortex liquid we may follow 
the procedure described above and find the temperature 
T cm (o-xiP) a t which C n '^ n in Eq. (|16f) for |An| = 1 is re- 
duced by a factor of 1/2 with respect to its zero temper- 
ature (r — > oo ) limit. This yields in the strong pinning 
limit, p 3> 1 (see Fig. 0J), the linear dependence 



T cm (a x ,p) ~ -86 x ApT a 



(19) 



Beside its influence on the vortex lattice melting tran- 
sition, the pinning potential can change the vortex lat- 
tice structure, both continuously and discontinuously. 
The zero temperature limit was discussed in Scc dlll All . 
Above the critical value p = 1 the lateral vortex positions 
coincide with the pin square lattice positions, £j = 0. For 
decreasing pinning strength below p = 1 , the configura- 
tion of the vortex lattice deviates continuously from the 
square structure to a lattice with vortices shifted along 
chains away from the pinning centers. 

Similar second order SPS to TPS phase transition (as a 
function of p ) is expected at finite temperatures. Indeed, 
as shown in Sec. ljlll A|) . the free energy functional /□ in 
Eq. @ is minimized at the stationary points = ^ c k = 
( — 1) 6>o , with cos 6*o = p for p < 1 , and at £ c fc = 
for p > 1. Thus, expanding /□ as a Taylor series in 
(£fc — £ c fc) about its stationary points it is clear that for 
p > 1 (when £ c fc = ) the expansion includes only even 
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powers ( due to the symmetry of /□ with respect to 
£k — ► — £j. ). Thus, at any finite temperature, T , the 
thermal mean values are equal to zero for p > 1 
, implying that for pinning strengths above the critical 
value p = 1 the mean vortex positions coincide with the 
square pin lattice. 



IV. PHASE DIAGRAM FOR THE LOWEST 
MATCHING FIELD 

The results of the previous sections enable us to draw 
a rather clear picture of the Vq — T -phase diagram, as 
shown in Fig. [5] In the strong pinning regime, p > 1, 
the pinning strength is so large that the gain in com- 
mensurate energy is larger than the vortex-vortex en- 
ergy gain at any temperature, and so the floating solid 
phase is not favorable. Thus, the vortex lattice melting 
in this region should take place directly from the SPS 
to the liquid phase, as described by the asymptotic ex- 
pression, Ea. l|19fl . which is equivalent to the straight line 
p w 0.29T/Tn in the large p regime of the phase diagram. 

In the small p regime the stable phase at low tempera- 
tures is the FS. Here the energy gain associated with cre- 
ation of the closed packed equilateral triangular vortex 
lattice, exceeds the energy cost of the incommensurate 
state. This state remains stable up to a relatively high 
temperature T ~ 2.8Tb , above which it melts into a 
vortex liquid state. The phase boundary in this region is 
vertical (i.e. independent of p ) since it is determined by 
the vortex- vortex coupling and not by the pinning energy 
(which is a constant in the floating state). 

In the low temperatures region of the phase diagram 
our analysis shows the existence of two phase transitions: 
At small pinning, increasing p above p c ks .25 transforms 
the FS discontinuously to a pin solid since the energy 
gain associated with the commensurate pin vortex solid 
exceeds the energy cost of distorting the closed packed 
equilateral triangular vortex lattice. The discontinuous 
nature of this transition is due to the fact that even in- 
finitesimal deviation from a commesurate configuration 
rises the pinning energy by a finite amount (i.e. at least 
from .6v to v). 

It turns out that the pin vortex crystal just above the 
commensurate-incommensurate transition is not a square 
lattice, as found by Reichhardt et alm&, but a triangu- 
lar one, with a unit cell which depends on the pinning 
strength. At T = it is a parallelogram with equal base 
and height, which transforms continuously to a square 
at p = 1. Similar continuous transition from a frus- 
trated triangular pin lattice to the SPS takes place at 
the critical pinning strength p = 1 at any temperature T. 
Interestingly, the corresponding horizontal transition line 
intersects the extrapolated SPS-L boundary line at T w 
3.44T n , p — 1 , that is in the close vicinity of the inter- 
section between the vertical FS melting line, T ~ 2.8T n , 
and the SPS-TPS line. 

It is not exactly known, however, how the FS-TPS 
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FIG. 5: Vo — T phase diagram: Solid lines - first order phase 
transitions from SPS to L phase (at large pinning strength), 
from TPS to FS (near zero temperature), and between FS 
and L phases. Dashed line - second order phase transition 
between partially pinned (TPS) and fully pinned (SPS) vortex 
crystals. The dashed-dotted line connects smoothly between 
the asymptotic SPS-L line and the low temperatures TPS-FS 
line (see text for explanation). 



boundary is extended beyond the zero temperature re- 
gion. It is conceivable that its high temperature sector 
coincides with the low temperature sector of the SPS 
melting line. This is due to the fact that, at a fixed 
value of p, the driving force for both transitions are ther- 
mal fluctuations involving sliding vortex chains, which 
suppress the pinning energy gained in the commensurate 
phase (i.e. the term — 4pcos(£fc) in Eq.©). In the SPS- 
L transitions, where the vortex- vortex interaction is rela- 
tively small, this suppression leads to uncorrelated vortex 
chains, resulting in melting. In the TPS-FS transitions, 
where the the vortex-vortex coupling is relatively large, 
the suppression of the pinning energy results in mutually 
correlated vortex chains, which lose correlation with the 
underlying pinning lattice. 

An intermediate pin solid phase of a triangular form 
has been also found in the London model calculation re- 
ported by Pogosov et alM*. However, in contrast to the 
Ginzburg-Landau model, discussed here, they predicted 
the vortex configurations shown in Fig. |3Jd as prefer- 
able below some critical value of the pinning potential 
strength. Above this value the symmetry of vortex lat- 
tice is changed discontinuously to the square symmetry 
of the pin lattice. 

Our proposed phase diagram, shown in Fig[3 thus con- 
sists of 2 pin solid phases, a floating solid and a liquid 
phase, dclimitted by 4 interphase boundary lines, which 
intersect at two nearby triple points. This result is sim- 
ilar to the phase diagram found by Reichhardt et al~, 
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using molecular dynamics simulations. However, the 
intermediate TPS phase obtained in our calculation, is 
missing in Rcichhardt et al. This seems to be due to the 
square boundary conditions imposed in the latter calcu- 
lation. Another difference concerns the zero temperature 
limit of the PS-FS line, which seems to approach p = 
in Reichhardt et al. 



V. CONCLUSIONS 

The influence of a periodic pinning potential on the 
vortex state of a 2D superconductor at temperatures 
well below the mean field T c has been studied within 
the framework of the GL functional integral approach. 
It is shown that attractive interaction of vortex cores 
to a commensurate pin lattice stabilizes vortex solid 
phases with long range positional order against violent 
shear fluctuations along the principal crystallographic 
axis. Exploiting a simple analytical approach we draw 
a rather detailed picture of the relevant vortex state 
p — T (pinning strength-temperature) phase diagram. 
In agreement with previous numerical simulations^, we 
have found a pinned, commensurate solid phase in the 
strong pinning-low temperature part of the phase di- 
agram, which melts into a vortex liquid at high tem- 
peratures, and transforms into a floating (incommensu- 



rate) solid at low temperatures. We have shown that 
at low temperature, similar to Ref^, there is an inter- 
mediate triangular phase, where vortices detaching from 
pinning centers remain strongly correlated with them. 
This pinned (frustrated) triangular solid transforms con- 
tinuously into the fully pinned (square) solid phase at 
p = I, and discontinuously to a floating solid at small 
pinning strengths. The zero temperature limit of this 
commensurate-incommensurate transition line occurs at 
a finite pinning strength (jp = p c m .25). 

It is predicted that for sufficiently clean samples, 
where random pinning is weak enough, application of 
an artificial periodic pinning array with an appropriate 
strength would stabilize a weakly pinned vortex solid 
phase with long range positional order. Exploiting the 
SANS method to the sample under these conditions 
one could therefore directly detect the shear stiffness 
anisotropy characterizing the ideal vortex lattice. 
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